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ABSTRACT 

Context. In order to accurately model giant planets, a whole set of observational constraints is needed. As the conventional constraints 
for extrasolar planets like mass, radius, and temperature allow for a large number of acceptable models, a new planetary parameter is 
desirable in order to further constrain planetary models. Such a parameter may be the tidal Love number k2. 
Aims. In this paper we aim to study the capability of to reveal further information about the interior structure of a planet. 
Methods. With theoretical planetary models we investigate how the tidal Love number k2 responds to the internal density distribution 
of a planet. In particular, we demonstrate the effect of the degeneracy of k2 due to a density discontinuity in the envelope of a three- 
layer planetary model. 

Results. The effect of a possible outer density discontinuity masks the effect of the core mass on the Love number ^2- Hence, there is 
no unique relationship between the Love number k2 and the core mass of a planet. We show that the degeneracy of k2 with respect to a 
layer boundary in the envelope also occurs in existing planets, e.g. Saturn and the Hot Neptune GJ436b. As a result of the degeneracy, 
the planetary parameter k2 cannot be used to further constrain Saturnian models and for GJ436b only a maximum possible core mass 
can be derived from a given k2. To significantly narrow the uncertainty about the core mass of GJ436b the combined knowledge of ^2 
and atmospheric metallicity and temperature profile is necessary. 
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1. Introduction 



Planets and satellites: individual: Saturn - Planets and satellites: individual: GJ436b - 



Ever sinc e the first extrasolar plan et around a solar-type star was 
detected dMavor & Ouelozlll995h . questions about the compo- 
sition and origin of extrasolar planetary objects (exo-planets) 
have been of major interest. Models of exoplanets are often lit- 
tle constrained based on the observable parameters mass, ra- 
dius, and effective temperature, in particular metal-rich planets 
d Adams et al.ll2008l) . For solar planets additional constraints are 
provided by the gravitational moments which have been mea- 
sured by spacecraft or Earth-based observations of the motion 
of satellites and hence are not accessible for extrasolar planets. 
However, a similar quantity does exist: the tidal Love number 
^2- To first order in the dimensionless number that describes the 
effect of rigid rotation or degr ee 2 tidal distortion k2 is equivalent 
to J2 (see e.g. lHubbardll984 . 

The tidal Love number kj, is a potentially observable param- 
eter. iRagozzine & WolU ( l2009h showed that the dominant source 
of apsidal precession of Hot Jupiters is the tidal interaction be- 
tween the planet and its star. This tidally induced apsidal preces- 
sion creates a unique variation in the transit light curve which is 
detectable by space-based missions like Kepler. Another possi- 
bility of determining k2 is the measurement of the orbital param- 
eters of a two-planet system in apsidal alignment ( iBatvgin et alj 
l2009l) . Due to tidal dissipation a coplanar two-planet system can 
evolve into a tidal fixed p oint which is ch aracterized by the align- 
ment of the apsidal lines (lMardlind2007h and both orbits precess 



* e-mail: ulrike.kraimn2@uni-rostock.de 



with the same rate. iBatv gin et al.l(l2009h showed that in this state 
the Love number k2 is a function of the inner planet's eccentric- 
ity. 

Like J2 for the solar system planets, k2, if known, can be 
used to further constrain the models of extrasolar planets as 
it is sensitive to the internal density distribution of the planet. 
Understanding the planetary interior is important for determin- 
ing not only physical processes but also the formation history. 
Hence, it is crucial to analyze what information can be extracted 
from a measured ^2 and its implications on the planetary interior. 

First, we will describe the definition and calculation of the 
Love numbers in Sect. [2] We also confirm the correlation be- 
tween the central condensation of a planet and its Love number 
k2 within a simple two-layer model. In Sect. [3] we introduce a 
more sophisticated three-layer planetary model and demonstrate 
the degeneracy of k2 with respect to the density discontinuity in 
the envelope. We apply these results to Saturn and to the Hot 
Neptune GJ 436b in Sect. |4] The main results of this paper are 
summarized in Sect. |5] 



2. The Love number k2 

2.1. Definition & calculation 

Love Numbers quantify the deformation of the gravity field of 
a planet in response to an external perturbing body of mass M, 
which can be the parent star, another planet or a satellite. M, 
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moving in a circular orbit of radius a around a p lanet, causes a 
tide-raising potential (IZharkov & Trubitsvnlll978h 



Wis) = 2 W„ = (GM/a)2(i/fl)"P„(cos0') 



(1) 



where s is the radial coordinate of the point under consideration 
inside the planet, ff the angle between the planetary mass ele- 
ment at s and the center of mass of M at a, and P„ are Legendre 
polynomials. Due to the tidally induced mass shift the planet's 
potential chan ges by Vi^'^i s) - Kn{s)W„{s), where K„(s) is the 
Love function dLovd 1 9 llh . Thus, at the planet's surface the def- 
inition of the Love numbers k„ reads 



Vf{R^) = kMR^) 



(2) 



As we are interested in low eccentricity synchronous or- 
bits, we concentrate on the purely hydrostatic tides. For the 
calculation of the Love numbers we follow the approach 
by [Zh arko v & Trubitsyn ( 1978), see also Gavrilov & Zharkovl 
(1 19771) and iGavrilov et all d 19751) . A Love number of degree n 
is obtained from 



Rpgo 



(3) 



where r„(/?p) is the value of the function T„{s) at the planet's 
surface, Rp is the radius of the planet and the surface grav- 
ity for the unperturbed planet. The function r„(i) satisfies the 
following second order differential equation: 



AnGp'is) n{n + 1) 
V'(s) ^2 



Tn{s) = 



(4) 



The radial coordinate is represented by s and p{s) and V{s) give 
the unperturbed density distribution and potential of the planet, 
respectively. The primes denote first and second differentiation 
with respect to the radius s. If the planet has an internal density 
jump then the jump condition for the function T is 



AnG 



T„{b^) = T„{b-) 
\p{b-)-p(b^)]T„{b) 



(5) 
(6) 



Here b is the radial position of the density jump and the "-n" 
("-") denotes a place just outside (inside) the discontinuity. This 
procedure assumes linear response, i.e. the tidal distortion is as- 
sumed to be small. Note that the only input needed for the calcu- 
lation of a Love number k„ is the radial density distribution p{s) 
of the planet. Hence, the Love numbers contain important infor- 
mation about the interior structure of a planet. In this paper we 
focus on the Love number ko (twice the apsidal motion constant 
described bv lSternel(ll939h ). 

2.2. Central condensation & two-layer model 

The Love number k2 is of special interest. It is a measure for 
the level of central condensation of an object: The more ho- 
mogeneous the planet in mass distribution, the bigger the Love 
number ko. Maximum homogeneity is represented by a planet of 
constant density, yielding the maximum value of ^2 of 1.5. If the 
planet is more centrally condensed, the Love number decreases. 
A planet with a density distribution of a « = 1 polytrope has a 
value of k2 = l5/n^ - 1 = 0.5198 . . . Planets with a core can 
have an even smaller Love number due to a stronger central con- 
densation, e.g. Saturn interior models with a ~ 10 core give 



k2 - 0.32. That is why it has been suggested to use k2 to infer 
the presence of a massive core dRag ozzine & Woill2009l) . 

First, we investigate the dependence of ^2 on the density dis- 
tribution within a simple two-layer model, consisting of a core 
with constant density A and an « = 1 polytropic envelope: 



Xr < X < 1 



sin^(l - x) 
qx 

P2(x) — A , < X < Xc 



(7) 
(8) 



The functions pi{x) and p2{x) give the radial density distribution 
of the envelope and core, respectively. The radius coordinate has 
been scaled such that x = 1 is the surface of the planet and 
Xc gives the radius of the core in planet radii (x - r/Rp). The 
core density A is always greater than pi(xc). The free parameters 
of this model are the core radius x^ and the quantity q. While 
Xc only influences the core density, q influences both the core 
density and the density in the envelope. The core density A is 
determined by the choice of Xc and q under the condition that 
just outside of the core hydrostatic equilibrium is satisfied. Due 
to the strong influence of q on the core density A, a change in q 
is equivalent with changing the ratio of core mass to total mass 
A/core /A/totaCl It is ucccssary to choose q < n since we must have 
pi(x) > 0. At the surface pi(l) = 0. 

The use of an « = 1 polytrope is motivated by two consid- 
erations. First, it is a good approximation for the behavior of 
hydrogen-helium mixtures for the range of pressures and tem- 
peratures of interest and is even still roughly correct if modest 
amounts of ice or rock are mixed with hydrogen and helium. 
Second, it yields an analytic form (given above) for the den- 
sity profile, thus making the Love number calculations straight- 
forward. The parameter q is directly related to the proportion- 
ality constant K in the assumed polytropic equation of state 
P = Kp^^^^^"\ Since our intent here is to understand the gen- 
eral nature of the dependence of the Love number on density 
structure and not to derive highly specific and precise models, it 
serves our purpose well. 

By choosing Xc and q we can construct different two-layer 
models. Figure [T] shows the Love number ^2 as function of 
A^core/Motai for diff'erent core radii of respectively Xc - 0.05, 
0.10, 0.15 and 0.20, and in the core-less case. To obtain different 
values for Mcoi-e/A^totai the parameter q has been varied. A value 
of ^ = TT gives the minimum ratio Mcore/A^totai as the core density 
A decreases with q. Reducing the q parameter causes a rapid in- 
crease of the core density A but only a modest increase of pi(xc). 
For the q parameter range of 0.97r < q < n the ratio of the en- 
velope density at the core-mantle boundary and the core density 
has values of pi(Xc)/A > 0.01. Models with too low q become 
unrealistic as they have an enormously small pi(Xc)/A. Note that 
the choice of the core radius Xc has a compensating effect and q 
could be choosen lower if Xc is bigger. In reality, lower k2 can 
be obtained in particular if the mantle density distribution does 
not follow an n = 1 polytrope and a compressive equation of 
state for the core material is used. As can be seen from Fig. [T] 
k2 decreases with increasing ratio of core mass to total mass. 
This behavior applies for all choosen core radii. For a planet of 
a given ratio of core mass to total mass, a bigger core radius Xc 
implies a lower density for that core and hence less central con- 
densation and a larger Love number, as Fig. [1] shows. The data 
points for the smallest mass ratios have been generated with the 
maximum value of q - tt. This means that for a two-layer model 
with a fixed core radius there is a minimum core mass. 



: 4nAxl/3 and Mtotai = A^corc + 4;r£ pi(x)x^ dx 
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Fig. 1. Love number k2 in dependence on the ratio of core mass 
to total mass for theoretical two-layer models consisting of a 
constant density core and an n = 1 polytropic envelope, see 
equations (|7]i and (jSj. The lines show solutions for core radii 
of jcc = 0.05 (solid), 0.10 (dashed), 0.15 (dotted) and 0.20 (dash- 
dotted). The black diamond shows the solution for the core-less 
polytropic planet. To obtain different values for Mcoie/Mtotai the 
parameter q has been varied, q = n gives the minimum value of 
A^cme/A^totai- Thick lincs show the region where Q.9n < q < n 
which corresponds to realistic pi(xc)/A ratios. 



These calculations within a simple two-layer model give an 
intuitive interpretation of the characteristic of k2 of being a mea- 
sure for the level of central condensation: the bigger the core 
mass, the more centrally condensed and hence the smaller the 
Love number kn. However, as we will show in the next section 
this simple deduction is no longer possible when there is another 
density discontinuity in the envelope of the planet. 

3. Degeneracy in three-layer models 

In this section we will investigate the behavior of k2 in a more 
complicated model. In addition to the density discontinuity at the 
core-mantle boundary we introduce another discontinuity in the 
envelope of the planet. A three layer structure is a common as- 
sumption in planet modeling and has been used f or mode l ing th e 
solar system giants Jupite r and Saturn (see e.g. iGuillotI (Il999l) . 
ISaumon & Guilio3 Jlool v Such a separation of layers in the 
planetary e nvelope can occur as a result o f demixing of hydrogen 
and helium jStevenson & SalpeteJlQ??!) . It could also arise from 
double diffusive convection in the presence of a density gradient 
that is in troduced during accretion or because of subsequent core 
erosion ( IStevensonll982h . Layer boundari es are also compatibl e 
with standard models of planet formation jHubbard et al.lll995l) . 
We define our theoretical three-layer model as follows: 



sin^i(l - x) 



Piix) = B 



qix 
singafa - x) 
qix 
P3(x) = A 



, < X < \ 

•)Cq JC 
< X < JCc ■ 



(9) 

(10) 
(11) 



It consists of a core with constant density A and two poly- 
tropic envelopes described by p\{x) (outer envelope) and piix) 
(inner envelope). The same characteristics as for the two-layer 
model apply. The layer boundary in the envelope is placed at x^. 




Xn 



Fig. 2. Lines of equal k2 (blue) and equal Mcore/A^totai (red) in 
Ap-Xm parameter space for a theoretical three-layer model, see 
equations (|9|, ( fTOb and ( fTTI ) with the fixed parameters qi - 
0.98;:, q2 = Q-99n and x^ =0.1. Numbers give the correspond- 
ing values of k2 and Mcoie/Mtotab respectively. The equi-A:2-lines 
demonstrate the degeneracy of k2 with respect to the position 
(x„,) and size (Ap) of a density discontinuity in the envelope and 
the ratio of core mass to total mass (a blue line intersects several 
red lines). 



We choose the parameters qi and q2 (near but smaller than n and 
not much different from each other). The location of the den- 
sity discontinuities x^ and Xm are also free parameters. The core 
density A is determined as for the two-layer model. The conti- 
nuity of pressure and gravity at Xm is used to calculate Xa. We 
can then choose B to get a specified non-dimensionalized size of 
the density jump Ap - (p2 - pi)/pi\^^^ . Summarizing, in this 
three-layer model we can vary the parameters Xc, x^, qi, q2 and 
Ap while A, B and Xa are determined by the parameters chosen. 

3. 1 . Love number and core mass 

The parameters characterizing the density discontinuity in the 
envelope are x^ and Ap, position and size of the discontinuity. 
Thus, in order to investigate the influence of the outer discon- 
tinuity on the Love number k2 we varied the parameters Xm and 
Ap from 0.5 to 0.9 planet radii (increment 0.01) and from 0.01 to 
0.5 (increment 0.01), respectively. For the example we give here, 
the other free parameters are fixed at = 0.987r, q2 - 0.99;r and 
Xc = 0.1. This choice is in order to keep the envelope structure 
close to an « = 1 polytrope. Together with a moderate value for 
the core size this example for the density distribution mimics a 
Jupiter-like planet. 

The change of the Love number k2 in dependence on the pa- 
rameters of the outer density discontinuity is shown by Fig. |2] 
Lines of equal k2 and their values are given. This demonstrates 
a degeneracy of the Love number ko with respect to the outer 
density discontinuity. One can always find many different x^,- 
Ap-pairs (that means different three-layer planetary models) that 
give the same k2. These models lie on one of the equi-A:2-lines 
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and hence one cannot distinguish between these models by a 
measurement of k2. In addition to the Love number k2 we cal- 
culated the ratio of core mass to total mass Mcoie/Motai in the 
same parameter space. Lines of equal mass ratio are also shown 
in Fig. |2] For a specific value of k2 planetary models with dif- 
ferent mass ratios are possible, see intersections of equi-A:2- 
and equi-Mcoie/Motai-lines. (Note that the resulting Mcoie/Motai 
have a very limited range because only and Ap change and 
A^core /A/total IS morc Sensitive to the q parameters which are fixed 
in our example here. Anyhow, other choices of the q parameters 
and/or would yield qualitatively similar results.) In the same 
manner, planetary models with identical mass ratios can have 
different k2 values. This demonstrates that a unique relationship 
between the core mass and the Love number k2 is no longer valid 
in a three-layer model. It can also be seen from Fig.|2]that plan- 
ets with a more massive core can have a bigger Love number 
than planets with a smaller core (compare e.g. left lower part 
with right upper part of Fig.|2]i. This is in contrast to the results 
that were suggested by the simple two-layer model where more 
massive cores lead to smaller k2. In conclusion, the effect of a 
possible density discontinuity in the envelope of a planet on the 
Love number ^2 can mask the effect of the core mass. 

Even though we found that there is a non-unique relation be- 
tween the Love number k2 and the core mass, Fig.|2]is still con- 
sistent with the assumption of k2 being a measure for the level 
of central condensation of a planet as the term 'central conden- 
sation' cannot be defined by the core mass alone. A big density 
jump Ap means that there is a significantly higher density in the 
inner than in the outer envelope. This constitutes a higher level 
of central condensation than an envelope without a discontinu- 
ity. That is why ^2 is generally smaller in regions with bigger Ap. 
On the other hand, a very small Ap means that there is a small 
difference in density at the layer boundarjQ. Hence, the planet 
is more homogeneous and the Love number ^2 tends to bigger 
values. In the region of small Ap the equi-^2 lines also become 
flatter, because the position x,,, of the discontinuity becomes less 
important for Ap — > 0, as the distinguishability between a two- 
layer and a three-layer model disappears. 




Fig. 3. Lines of equal moment of inertia Cnd for the same three- 
layer model as in 93. H and in the same Ap-Xm parameter space. 
Blue lines are for the moments of inertia obtained directly from 
the Love numbers k2 via equation ( fTSl ). Green lines show the 
moments of inertia calculated from the density distribution with 
equation ( fT4l i. 



Since the tidal potential ([T]) has the same functional form 
as the rotational potential, the tidal response of a planet can be 
treated analogously to the rotational one. For a liquid or perfectly 
elastic planet the tidal response coefficient Aj^ = k2/3 is the 

same as the rotational one (lHubbardll984 . Substituting A^^ for 

. (1) 



Aj Q in equation (fT2l l. we find: 



3.2. Moment of inertia 

Another quantity closely related to the structure of a planet is 
its moment of inertia C, or its non-dimensional form Cnd - 
C/MRp. While it has been possible to obtain moments of in- 
ertia fromjneasurements of precession for Earth and Mars (see 
e.g. iFoUcner et alj|1997h . it is more difficult for the giant plan- 
ets. For the modeling of the internal structure of the giants, the 
gravitational moments have been used instead. The first-order 
response of a planet to rotational distortion (quantified by ^2) 
can be related to the planet's nondimen sional polar mo ment of 
inertia with the Radau-Darwin equation (lHubbardll984t) : 



C 



(RD) 
nd 



3A^'''o + 1 



- 1 



(12) 



where M and Rp are the mass and the radius of the planet, re- 



spectively, and to linear order A*'g - J2/q with q - oj'^Rp/GM, 
the ratio of centrifugal force to gravity. This equation implies 
that the response coefficient Aj g is a function of the moment of 
inertia only. 



^ The limiting case Ap 
was discussed above. 



is equal with a two-layer model like it 



c. 



(RD) 



SXko + l 



- 1 



(13) 



With this relation we can calculate the nondimensional mo- 
ment of inertia from our previously generated values of the Love 
number ^2- Another, more direct, method to obtain the moment 
of inertia for our models is to integrate over the density distribu- 
tion: 



I Jq x*p{x) dx 
j^' x^p(x) Ax 



(14) 



We calculated moments of inertia for the reference three- 
layer model of 93.11 In analogy to Fig. [2] the characteristics of 
the envelope layer boundary are varied. Fig. [3] shows lines of 
equal moment of inertia, obtained from equation ( fT3] l or (fT4l) . 

Of course the equi-Cnd-lines generated from ( flj] ) are equiv- 
alent to the equi-fc2-lines of Fig. |2] as (fT3b is a one-to-one cor- 
respondence. Obviously, they are different from the results ob- 
tained from (fT4l) . This is due to the approximations made for the 
derivation of equation ( fT3] ). 

The Radau-Darwin equation (fT2l) itself is an approximation. 
It only approaches a correct result Cnd — > 0.4 for the limit of 
a completeley homogeneous planet. It is not true in the limit 
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Cnd (see e.g. lMurrav & Dermotll[T999l) . This feature of the 
Radau-Darwin equation can be seen in Fig. [3] The difference 
between the equi-Cnd-lines becomes bigger with increasing Ap, 
corresponding to stronger central condensation. 

We conclude that the Radau-Darwin relation is a too crude 
approximation to describe the moment of inertia of gas giant 
planets. 



4. Application to existing planets 

In this section we apply the results obtained for the simple mod- 
els shown above to interprete interior models for real planets. 
We will show that the degeneracy of k2 is not just an artifact of 
our theoretical three-layer model but actually occurs in existing 
solar as well as extrasolar planets. 



4.1. Saturn 

As the second largest planet in our solar system, Saturn has 
attracted much attention in the past. With the help of space 
missions like e.g. Cassini observational data have been de- 
rived that provides important constraints for modelling Saturn. 
Voyager measurements have yielded a depletion of heli um com- 
pared to the solar abundance dConrath & Gautie^l2000l) . In con- 
clusion, the missing helium must be hidden deep inside the 
planet. The helium depletion in the outer part of the envelope 
can b e explained by a pha se se paration of hydrogen and he- 
lium. iLorenzen et al.l (l2009h and [Morales et al.l (|2009|) showed 
recently that demixing occurs on large scale at standard Saturn 
interior model conditions, motivating a layered-envelope as- 
sumption for Saturnian models. Furthermore, the demixing the- 
ory allows to reproduce the correct age of Saturn whic h was not 
possible without demixing (iFortnev & Hubb^dll2003h . 

As Saturn can be modeled in a three- or more layer approach, 
it is also affected by the A;2 -degeneracy. We dem onstrate this 
effect by examining selected Saturn models from iNettelmannI 
(l2009h . Those models have a rocky core based on the rock equa- 
tion of state (EOS) by Hubbard & Marley ( 1989) and two isen- 
tropic envelopes of hydrogen, helium and water (as represen- 
tative for metals) based o n the linear r rdxing Rostock equation 
of stat e ( LM-REOS), see iHolst et al.l (l2008l) . iKietzmann etaP 
(l2007h and lFrench et al.l (l2009b . respectivelv. The pressure at the 
layer boundary in the envelope is called the transition pressure 
f 12. By changing P[2 the layer boundary can be shifted. 

The variation of P12 affects the core mass of Saturn. We as- 
sumed a layer boundary between 1 and 3.8 Mbar, where P12 - 
1 Mbar gives a rocky core mass of 14 Me and P12 = 3.8 Mbar 
gives the maximum value of P12 as the core vanish es when the 
layer boundary is put that deep inside the planet dNettelmannl 
l2009l) . It is important to note that all these models reproduce the 
observational constraints, including the gravitational moments 
J2, Ja and J(,. This means that the parameter P\2 is not fixed. 
As a consequence, there is a large uncertainty in the core mass 
of at least 14 M© . We calculated the Love number k2 for these 
models and found for all these Saturn models ^2 - 0.32. Even 
though the core mass varies fom to 14 Me , no significant dif- 
ference in the k2 values can be seen. This is an expected result 
because k2 is proportional to J2 to first order in the expansion 
of the planet's potential (see e.g. iHubbardll 1984) and all models 
match the observed J2. 

The constant A:2-value demonstrates the effect described in 
the previous section: The influence of the outer density disconti- 
nuity masks the effect of the core mass. For one k2 value different 



planetary models can be found. Increasing the transition pressure 
P\2 is equivalent to shifting the density discontinuity inwards to 
smaller (compare theoretical three-layer model in sectionO. 
At the same time the size of the density jump Ap grows in order 
to obtain models that are consistent with the measured gravita- 
tional moments. So all Saturn models presented in this work lie 
on an equi-A:2-line (compare Fig.|2]i, showing the degeneracy of 
the Love number k2 with respect to the outer layer boundary in 
Saturn. 

Figure |2] also predicts a growing core mass with a deeper 
layer boundary and increasing density jump. However, this is 
not applicable here because our Saturn models are based on the 
additional constraints of matching the gravitational moments. 
These constraints were not considered for the simple models in 
Sect. |3] In the planet modelling procedure the gravitational mo- 
ments are adjusted by varying the envelope metallicities. This 
means that the metallicities in the envelope are determined by 
the observed gravitational moments for a specified P12. Shifting 
the layer boundary inwards (bigger P12, smaller x^) forces the 
metallicity of the inner envelope to increase in order to ensure 
consistency with the gravitational moments. These "additional" 
metals for the inner envelope are taken from the core, which 
is why the core mass shrinks with increasing transition pres- 
sure P12. 

4.2. G J 436b 

For modelling extrasolar planets the Love number ^2 can be a 
valuable constraint because the gravitational moment J2 is not 
accessib le to observatio ns. To first order ^2 is equivalent to J2 
(see e.g. lHubbardll984 and is potenti ally observable due to ap- 
sidal precession of the planet's orbit jRagozzine & Wol3l2009l) 
or measurement of the orbital parameters of a two-planet sys- 
tem in apsidal alignment (Mardling 2007). However, interior 
models of extrasolar planets are also subject to A:2 -degeneracy. 
Here we demonstrate this effect with the example of the Hot 
Neptune GJ 436b. Numerous investigations about this transiting 
planet have been made in the past (see e.g. lAdams et al. I l200i 
Figueira et al. 2009; Torres et al. 2008]). Based on mass, radius, 
and temperature measurements alone, a large variety of models 
is possible, ranging from a "water world" to a Super-Earth. We 
have performed exten sive calculations for mode lling the interior 
structure of GJ436b ( iNettelmann et al.ll20T0bl) . We considered 
two-layer models with a rocky core and one envelope where wa- 
ter is homogeneously mixed into H/He as well as three-layer 
models consisting of a rocky core, a water layer and an outer 
H/He layer with metal (water) abundance Zi . For selected mod- 
els we calculated t he Love number k 2, for details about the mod- 
elling see .Nettelmann et al. l (l20T0bl). Here we present a new and 
extended version of Fig. 2 in INettelmann et al. I dloTobI) in order 
to show the degeneracy of k2. 

Figure |4] shows Love numbers k2 and core masses of vari- 
ous two- and three-layer models of GJ 436b. All models shown 
here have an isothermal outer layer from 1 to 100 bar. For com- 
parison, Fig.|4]also contains the solar planets Jupiter, Saturn and 
Neptune and a 20 Me water planet. For two-layer models solu- 
tions are shown which have a metal abundance of Zi = 0, 0.5, 
and the largest possible value when Mcore = 0. Model series with 
three layers are displayed for constant outer envelope metallici- 
ties of Zi = 0, 0.5, and 0.7. Solutions are shown for two different 
surface temperatures of 700 and 1300K. Even though possible, 
we consider colder mod els to be unrealistic as recent results by 
ISpiegel et al.1 (1201 Oh and lLewis etaTI (1201 Oh predict 1300 K and 
1 100 K, respectively. In such warm models water is in the plasma 
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Fig. 4. Love number k2 and core masses of two-layer models 
(squares) and three-layer models (diamonds) of GJ 436b for tem- 
peratures of 700 K (cyan) and 1 300 K (orange) of the isothermal 
(1-lOObar) atmosphere. Color coded circles denote models of 
Jupiter (red), Saturn (yellow), Neptune (blue) and a 20 Me water 
planet (black). Three-layer models with the same metal abun- 
dance in the outer envelope are connected with dashed (Zi = 0), 
dotted (Zi = 0.5) and dash-dotted (Zi - 0.7) lines. The grey 
dashed line marks a highly degenerate area constructed to con- 
tain all Zi - Q three-layer models (the position of this line h as 
been shifted compared to Fig. 2 in iNettelmann et all ( 1201 Obi) in 
order to account for the new 1 300 K model series). 

phase and miscible with H/He. Hence, the water EOS applied to 
the inner envelope of such warm three-layer models represents a 
mean density of a H/He/H20/rock mixture. 

For GJ 436b we find 1, = - 0.72. Models of GJ 436b 
with k2 - 0.02 - 0.24 can have any core mass between and 
~ 0.9 planet masses. This means that a measurement of k2 in 
this regime would not help constraining the interior further. This 
highly degenerate area includes all Zi = three-layer models 
(note that this area has been extended compared to Fig. 2 in 
INettelmann et alJ (1201 Obh in order to account for the new 1 300 K 
model series). However, one can predict maximum possible core 
masses for specified k2. The uncertainty in the core mass de- 
creases with increasing k2. The upper limits of possible core 
masses, if k2 and the atmospheric temperature profile are known, 
are necessarily given by two-layer models with one homoge- 
neous envelope. Any redistribution of metals from an outer to an 
inner part of the envelope would, if the central condensation is 
to be kept at constant k2 value, require a compensating transport 
of core material to the envelope. The core mass would therefore 
decrease and drop below these upper limits. This is shown by the 
three-layer models with constant envelope metallicity of respec- 
tively Zi = 0, 0.5, 0.7, which all have lower core masses than the 
two-layer models. Three-layer models contain a density discon- 
tinuity in the envelope, dividing it into an inner water envelope 
(Z2 = 1) and an outer H/He/H20 envelope (Zi = 0, 0.5 and 0.7). 
The three-layer model approach gives the additional free param- 
eter P12, which is the pressure at the envelope layer boundary. 
To construct the equal-metallicity lines the metallicities are kept 
at constant values (Z2 - 1 and Zi - 0, 0.5, 0.7) and f 12 is 
varied from a minimum value where M^ore = to a maximum 
value where the water layer vanishes and the model is identical 
to a two-layer model with the same envelope rnetallicity (for the 
modelling procedure see also INettelmann et al.l (l2010bl) V 



Models of equal envelope metallicity Z\ show systematic 
behavior in M^oK-k2 space. For low outer envelope metallici- 
ties the Love Number k2 shrinks with growing core mass due 
to the stronger central condensation created by the bigger core. 
With a high metal abundance in the outer envelope this behavior 
changes and k2 increases slightly with growing core mass. The 
higher amount of metals in the outer envelope leads to a smaller 
density gap at the envelope layer boundary so that the planet be- 
comes more homogeneous than planets with a low metallicity 
envelope. As can be seen from Fig. HI this effect, and especially 
the transition between the two different slopes, is strongly tem- 
perature dependent. For lower temperatures (see 700 K curves) 
k2 increases with the core mass already at smaller metallicities 
than for warmer temperatures (compare 1300K curves). That is 
because colder envelope temperatures increase their density and 
hence the homogeneity of the planet. The upper limit of possi- 
ble core masses increases with decreasing atmospheric tempera- 
tures. 

Concluding, an observational k2 would imply a maximum 
possible core mass. But, especially for ^2 < 0.24, it would 
not help to further constrain interior models of GJ436b be- 
cause in that regime the solutions are too degenerate. However, 
ak2 > 0.24, if measured, would indicate a maximum core mass 
Mc < 0.5 Mp and large outer envelope metallicities. Only with 
additional information about the planet's atmospheric compo- 
sition (Zi) together with ^2 interior models of GJ436b can be 
further constrained. If we had knowledge about the atmospheric 
temperature profile and metallicity and the Love number k2, the 
core mass could be determined or at least strongly narrowed, de- 
pending on the slope of the equi-metallicity-lines. However, not 
all aspects of metallicity are potentially measurable in an exo- 
planet, only those that are the consequence of elements which 
are in the gas phase in the observable part of the atmosphere. As 
an example in our own solar system, the metallicities of Jupiter 
and Saturn are not measurable spectroscopically because water 
forms clouds deep below the photosphere. 

The results for GJ 436b are in perfect agreement with our re- 
sults obtained from Fig. [2] For a given k2 value several planetary 
models and core masses are possible (highly degenerate area in 
Fig. lU. Even for constant outer envelope metallicity there are 
also models that have a larger k2, indicating a rather homoge- 
neous planet, but nevertheless a more massive core. This behav- 
ior is caused by the outer density discontinuity at the water-H/He 
layer boundary as discussed in Sect.[3]and predicted by Fig.|2l 

5. Summary 

In this paper we investigated the eff^ect of the density distribu- 
tion of a planet on its tidal Love number k2 in order to find out 
what conclusions can be drawn from a measured k2 on the in- 
ternal structure of a planet. We confirmed that the Love number 
k2 is a measure of the level of central condensation of a planet. 
However, in a three layer model approach ^2 is not a unique func- 
tion of the core mass. There is a degeneracy of ^2 with respect to 
a density discontinuity in the envelope of the planet. It is possi- 
ble to have several acceptable models for a given ^2 value, which 
can differ significantly in core mass. The effect of the outer den- 
sity discontinuity on ^2 is compensating the effect of the core. 
Furthermore, we showed that the Radau-Darwin relation is a too 
crude approxomation to describe the moment of inertia of gas 
giant planets. 

We verified our results on k2 with models of existing planets. 
For Saturn the freedom to place the layer boundary in the enve- 
lope leads to a high uncertainty in the core mass. Regardless of 
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the core mass all Saturn models have the same ki, demonstrating 
the degeneracy caused by the outer layer boundary. 

For extrasolar planets the Love number ^2 can be an equiv- 
alent constraint to J2 for the solar system planets. However, one 
has to be careful with estimates about the core mass derived 
from k2 as degeneracy may also occur in extrasolar planets. For 
GJ436b we find a highly degenerate area of k2 < 0.24 where 
a measurement of would barely help to further constrain the 
interior models. Only a maximum possible core mass and for 
k2 > 0.24 a large metallicity can be inferred. With additional 
knowledge about the atmospheric metal abundance the uncer- 
tainty about the core mass could be significantly narrowed. 

Tabulating ^2 values of various planetary models can 
prove to be very useful once ^2 is actually measured for ex- 
trasolar transiting pla nets. For ins tance, for the Super-Earth 
GJ 1214b Nettelmann et all (1201 Oah demonstrated that H/He or 
water envelopes result in significantly different values of k2. 
Furthermore, we have shown in this paper that even though the 
Love number ^2 is a degenerate quantity it can help constraining 
the core mass of a planet. Knowledge about the core masses of 
planets is highly desired because it is thought to help to distin- 
guish between the possible planet form ation scenarios of core 
accretion (see e.g. Pollack et al.lll996l) and gravitational insta- 
bility (see e.g. iBossI 1997 ). However, one has to keep in mind 
that core accretion models can also result in very small cores 
of 1.7 Me or 0.25 Mg in the case of grain-free o r even metal- 
free envelopes, respectively dHori & Ikomall20Tol) . On the other 
hand, gravitational instability models allow the formation of a 
massive core as well if the protoplanet is cold en ough for grain 
settling to take place ( iHelled & Schubert! llOOSh . A clear dis- 
tinction between the two formation models can only be made 
for ma ssive extra-solar giant planets > 5 MT. lHelled & Schubert! 
(l2008h showed that such massive protoplanets formed by disk 
instability cannot build up a core at all due to their high internal 
temperatures and evaporation of the grains. 
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